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Abstract 

We study the B s —»• n + n^ and Bd —> K + K~ decays in the standard model and the family non-universal Z' 
model. Since none of the quarks in final states is the same as the initial quark, these decay modes can occur only 
via power-suppressed annihilation diagrams. Despite the consistence of the standard model prediction with the 
available data, there is a surviving room for a light Z' boson. Taking into account the Z' contribution, we find 
theoretical results for branching fractions can better accommodate the data. With the relevant data, we also 
derive a constraint on the parameter space for the Z'. Moreover, for the Bd —¥ K + K~, both the direct and the 
mixing-induced CP asymmetry are sensitive to the couplings between Z' and fermions in the parameter spaces 
constrained by data. The measurements at future experimental facilities, including the LHC-b, Bclle-II and the 
proposed high energy e + e _ collider, will provide us useful hints for direct searching for the light Z' boson. 


1 Introduction 

Since the discovery of the Higgs boson at the Large Hadron Collider (LHC) [Tj, the search for new physics (NP) 
degrees of freedom beyond Standard Model (SM) becomes one of the most important tasks of high energy particle 
physics. In many NP models, an extra U(l) r gauge symmetry is often introduced based on various motivations in 
new physics beyond SM, resulting in an additional massive neutral gauge boson usually called the Z' boson. Quite 
a few models are of this type, such as grand unified theories based on the gauge groups SO(10) [2], Eq model [3], 
supersymmetric models [3], and string inspired models [5] (for a review, see Ref. 0)- Although the 17(1)' charges 
are usually family-universal, it is not mandatory to be so, and the family non-universal Z' has been introduced in 
some models, such as in aforementioned Eq model [5]. 
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On the experiment side, many efforts have been expended to search for the Z' directly at the LEP, Tevatron, 
and LHC. With the assumption that the Z' couplings to the SM fermions are similar to those of the SM Z boson, 
the direct searches for the Z' can be performed in the dilepton events. At this stage, the lower mass limit has 
been set as 2.86 TeV at the 95% confidence level (CL) from collisions at 8 TeV with an integrated luminosity of 
19.5RC 1 by using e + e _ and p + p~ [7] events, and this value becomes 1.90 TeV using the t + t~ events [8|. 

However, if the Z' boson does not couple to leptons, the above constraint from the LHC is no longer valid. 
Theoretically, such leptophobic Z' boson can be realized in Eq model [3] , the phenomenological studies at the LHC 
has been recently explored in Ref. 0 ■ In complementary to the direct search, some characters of the leptophobic Z' 
boson can also be constrained indirectly from the “low” energy flavor physics. The family non-universal Z' boson 
may induce tree-level flavor changing neutral currents (FCNC) and thus they are severely bounded by experiment, 
most notably meson mixing [ID]. Other effects of the FCNC in flavor physics have also been studied in past decades 
mi mm- Motivated by the above arguments, we aim to in this work perform a comprehensive analysis of the 
impact of a family non-universal Z' boson on the pure annihilation decays B d K + K~ and B s —> . Since 

these modes are power suppressed in the heavy quark limit, their branching ratios are expected to be very small, 
and the sensitivity to NP can be then enhanced. 

Experimentally, the decay mode B s —> n + n~ was firstly reported by the CDF collaboration 

B(B S —>• 7r + 7r - ) = (0.57 ±0.15 ±0.10) x 10“ 6 [BO, (1) 

and it was soon confirmed by the LHCb collaboration with 0.37 fb _1 data as 

B(B S 7r + 7T“) = (0.95±g;f^ ± 0.13) x 10" 6 Qjg. (2) 

So, the averaged result is given as flGl : 

B(B S -> 7T+7T-) = (0.73 ± 0.14) x 10 -6 . (3) 

The branching fraction of another pure annihilation decay mode Bd —> K + K~ has been also measured as [161 : 

B{B d -+ K+K~) = (0.12 ± 0.06) x 10" 6 (4) 

Theoretically, within QCD factorization (QCDF) approach ITT] , only an order of magnitude estimate can be 
given for these two decays through introducing new phenomenological parameters (pa and 4>a) or an effective 
gluon propagator [18] due to the existence of the endpoint singularity. The predicted branching fractions are at 
the order of 10 -8 [JTJQj)]. Moreover, the effects of SU(3) asymmetry breaking have also been discussed in [20 1. 
On the contradiction, the perturbative QCD (PQCD) approach [TT] retains the transverse momenta of all inner 
quarks, and thus the endpoint singularity disappear. This makes the perturbative calculations of pure annihilation 
decay modes reliable. On the basis of PQCD, the decays B s —> and B d —> K + K~ have been explored in 

Refs. [2H123] and [2T], respectively. In Ref. j55J, the authors have revisited these two decays with new parameters 
(especially for the distribution amplitudes of light mesons), and the obtained results are in agreement with the 
experimental data well. Despite the agreement, by comparing the predictions of [25] with the experimental result, 
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one can find that the LHCb measurement has a central value larger than the theoretical result, which may indicate 
some room left for survival of a light Z' boson. In the following we will use the PQCD approach and investigate 
the impact of the family non-universal leptophobic Z' model on the CP asymmetries of these two decays. Our 
results can be stringently tested at the LHCb experiment, Belle-II, and future high energy e + e~ collider. 

This paper is organized as follows. In SecGl after a brief introduction to the PQCD approach, we will present 
the numerical results of B s —> and Bd K + K~ in SM. In Sec[51 we will discuss the effects of the Z' on the 

branching fractions and CP asymmetries of these two decay modes. At last, the conclusion will be drawn in the 
Sec0J 


2 SM Calculation 


In this section, we will start with the effective weak Hamiltonian for the b —> D (D = d, s) transitions, which are 
given by [25] 


Keff = 


Gp 

V2 


10 


E 

q=u,c 


V qb V* D (CiO? + c 2 o 2) - v tb v t * D J2c l o l \ + H.c., 


2—3 


(5) 


where V qb (D) are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The explicit expressions of the local 
four-quark operators Oi (i = 1,..., 10) and the corresponding wilson coefficients Ci at different scales have been 
given in Ref. [25]. Note that 0\ 2 are tree operators and others O3_io are penguin ones. 

The PQCD approach is based on the kx factorization, and has been applied to calculate the non-lcptonic B 
meson decays for many years mm- In this approach, the decay amplitude is conceptually written as 

A ~ f dxidx2dx3 I bidbib2db2b 3 db 3 Tr^C(t)$B(xi,bi)$2(x2,b2)$3(x3,b3)H(xi,bi,t)St(xi)e~ S (t) j, (6) 


where Xi are the momentum fractions taken by light quarks in each mesons, and bi are the conjugate variables 
of the transverse momenta of light quarks. “Tr” means the trace over both Dirac and color indices. In light of 
the factorization hypothesis, the wilson coefficient C{t) encapsulates the dynamics from mw down to the scale t, 
where t ~ 0(Mb/ 2) is the typical scale of the concerned annihilation type decays. The hard part H , involving the 
four-quark operators and the hard gluon, describes the hard dynamics characterized by the scale t, and it can be 
calculated perturbatively. The wave function standing for hadronization of the quark and anti-quark into the 
meson iff, is independent of the specific processes and thus universal. The factor St(xi) arises from the resummation 
of the large double logarithms (In 2 Xi) on the longitudinal direction, while the Sudakov form factor e~ s ^ is from 
the resummation of the double logarithm In 2 kx- Fortunately, the endpoint could be smeared effectively with the 
help of these two functions, which makes our calculation reliable. 

In particular, the wave functions («, /3 being Dirac indices) are decomposed in terms of the spin structure, 

la/3, 7 £a, (75<r A “ y ) a /3, (7 M 75)a/3 and (75)0 ,p. For the heavy pseudo-scalar meson B q (q = d,s ) meson, the wave 
function &B,a@ is given by 


®B, a p(x,b) 


{(/is75)a/3 + 75a/3)} <Mx, b), 


(7) 
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Figure 1: The Feynman diagrams for annihilation contribution, with possible four-quark operator insertions 


where N c = 3 is the color degree of freedom, and p B is the momentum of B meson. The scalar distribution 
amplitude (j) B is normalized by its own decay constant f B 

Ib 


/ (f> B (x, b = 0 )dx = 

Jo 


2 v / 2 W 


( 8 ) 


In this work, we employ the function 


4> B {x,b) = N b x 2 (1 


\2 

x) exp 


[ if 

^ xnriB N 

\ 2 ^ 1 

2 

V <XB , 

/ 2 


(9) 


where the shape parameter u> Bd = 0.4 GeV (uj Bs = 0.45 GeV) has been adopted in all previous analysis of exclusive 
B d ( s } meson decays J5Tj 1321 122] ■ 

In contrast to the heavy meson, the wave functions of light meson </>m are much complicated due to the non- 
negligible chiral mass. Taking the K + meson as an example for illustration, we define its wave function as 


®K+,ap(x, b) 


i 

V5K 


0 lbi>K)(!>K{x , b) + TO 0 if75^(a;, b) + vhqk^H - 1)<^ (x, b) 



( 10 ) 


where p B is its momentum, and = Bn\/(m u + m s ) is the aforementioned chiral mass, v and n are unit vectors, 
and v (n) is (anti-)parallcl to px- As nonperturbative parameters, the light cone distribution amplitudes (LCDAs) 
should be fixed by experimental data in principle. Though there is no direct experimental measurement for 
the moments yet up to now, the non-leptonic charmless B q decays already give much hints on them [sum. Since 
the PQCD approach had already given very good results for these decays, especially for the direct CP asymmetries 
in B° —y and B° —> K + ir~ decays, we shall adopt the well constrained LCDAs of the mesons in these 
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papers [27] (see [21] for a summary and update of the LCDAs): 


<t>ni x ) = 

^x(l-x)[l+0A4C 3/2 (t)}, 

<t> P (x) = 

^=[l + 0.43C 2 1/2 (f)], 

<t>Z( x ) = 

-^ [ C'i 1/2 (i) + 0-55C 3 1/2( t)], 

<(40 ) = 

^*(1 ~ *)[1 + 0.17<7 3/2 (:f) + 0.2Cf 2 (i)], 

<Pk(x) = 

^=[l + 0.24<2 /2 (i)], 

4>k( x ) = 

-^Pi 1/2 (i) + 0-35C 3 1/2( t)], 


with Gegenbauer polynomials defined as: 

Cl /2 (t) = t,Cf 2 (f) = 3f,C 2 1/2 (f) = i(3t 2 - l),C 3/2 (f) = ^(5f 2 - 1 ),C 3 1/2 (t) = if(5t 2 - 3), (12) 

1/2 3/2 

and t = 2a: — 1. It should be stressed that we have dropped the terms proportional to C 4 , and only kept the 
first two terms, following the arguments of [23] . 

Now we turn to calculate the hard part H. According to the effective Hamiltonian, eq.©, we can draw four 
kinds of Feynman diagrams contributing to the Bd —> K + K~ and B s —> decays at the leading order, as is 

shown in FigJT] The four diagrams are classed into two types: (a) and (6) are factorizable diagrams , and (c) and 
(■ d ) are non-factorizable ones. Due to the current conservation, the contributions from the factorizable diagrams 
(a) and (6) will be canceled exactly by each other, so that contributions from diagrams (a) and (b) are null. As 
concerned as diagrams (c) and (d), by inserting the possible operators, we can obtain the amplitudes for the non- 
factorizable annihilation diagram M^ r [ J n and M^ n , where LL stands for the contribution from (V — A)(V — A) 
operators, and SP for the contribution from (S — P){S + P) operators which result from the Fierz transformation 
of the (V — A) (V + A) operators. The expressions of M^ n and M^ n and the inner functions can be found in [23] . 
Finally, we obtain total decay amplitudes for concerned decays as 

A( B S 7r+TT-) = ^{v ub V u * s M^ n [C 2 ] - V tb V* (m“ [C 4 + C 10 ] + M a Z [C 6 + C 8 ] 

^ <->•*■+ + ^ann[C& ~ Ck]*- <-j.,r+) ) J (13) 

A(B d K+K-) = ^{v ub V; d M^ n [C 2 ] - V tb V t * d (M^ n [C 4 + C w ] + M s J n [C 6 + C s ] 

+ M ann\C 4 ~ -jCw ]k~++K+ + M ann[C 6 ~ 7} Cs]k-++K+^ } • (14) 

In ea. (Thill . when 7r + and n~ exchanging, M^ n obtain the results because of SU(2) symmetry. Furthermore, if we 
ignore the small x 4 (the momentum fraction of s quark in the B s meson) in the denominators, M^ n is same as 
47 0 S ]m,t° 0 Q. However, for Bd K + K~, Mann SP ^ do not share same formulaes with Mann SP ' > \K~<^K+ due to the 
difference between the mass of up (down) quark and that of the strange quark, and such difference might affect 
x In Ref. [25], there are typos in eqs.(27) and (28). 
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the direct CP asymmetry. In fact, in our calculations there are many uncertainties, the most important one of 
which is from the distribution amplitude of initial heavy meson, because it cannot calculated directly from QCD 
till now yet. In the following work, we shall vary the shape parameter u>B d = 0.40 ± 0.05 and u>b 3 = 0.50 ± 0.05. 

Furthermore, the contributions from next leading order (NLO) have not been done. In the current work, to estimate 

the uncertainties of NLO , we simply vary t from 0.8f to 1.2 1 , where t is the largest scale in each diagram and 
the expressions of them have been given in |Z5] . Combining all above uncertainties, we obtain the CP-averaged 
branching fractions of two decay modes 

B(B s ^7r + 7r-) = (5.5lo.c>) x 10 -7 . (15) 

B(B d ->K + K~) = (1.9+°;|) x 10~ 7 . (16) 


Since the uncertainties from the 7 r, K meson distribution amplitudes are very small, we will not discuss them here. 
In discussing the B meson decays, we usually define direct CP asymmetry as 

\A(B q ^f)\*-\A(B q ^7)\ 2 


A dir _ 1 ^ v - *— ^ 

= ■ — 


(17) 


\A{B q ^f)\* + \A(B q ^f)\* 

Moveover, because the final states tt + tt~ , K + K^ have definite CP-parity, one can measure the time-dependent 
decay width of the the B q —► / decay [29]: 


T(B(t) —> f) oc cosh + Hf sinh — Aq p cos(A mt) — Sf sin(Amt), 


ATt\ 


(18) 


where A m = mn — mi > 0 is the mass difference, and Ar = F# — is the difference of decay widths for 
the heavier and lighter P° mass eigenstates. Correspondingly, the time dependent decay width F (B° q (t) —> f) is 
obtained from the above expression by flipping the signs of the cos(A mt) and sin(Amf) terms. The Sf and Hf 
that can be extracted from the-time dependent decay width are defined as 

2Pe[A] 


_ 2/m [A] 

s ' = TTw 


with 


A = 


ri f e 


Hf ~ 1 + |A| 2 ’ 

2 i( A{B q f ) 

A(B q -+ f) ’ 


(19) 


( 20 ) 


where 77 / is +1(—1) for a CP-even (CP-odd) final state / and e = arg [—V c bVtqV* q V t * b }. In SM, the predicted results 
are listed as 

A d S r P {B s -S- 7T+7T-) = (-1.5 ± 0.2)%, 

Sf(B s ->• 7T+7T - ) = 0.11 ± 0.01, (21) 

H f (B s ^tt+tt") =0.99; 

A d S r P {B d -a K+K~) = (37^)%, 

S f {B d -+ K+K~) = -0.81 ±0.05, (22) 

U f (B d K+K~) = -0.45 ± 0.05, 

For B s —> , both branching fraction and CP asymmetry parameters agree with previous studies [2211231 [25] , 

and small differences are from the uncertainties of the CKM matrix elements. For B d —> K + K~ , our branching 
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fraction consist with prediction of [25| . but the direct CP asymmetry is much smaller than theirs because they 
may omitted the effect of the SU(3) asymmetry in the LCDAs of K meson. Compared with the experimental data, 
although our branching fractions are consistent with data after considering the uncertainties of both theoretical 
and experimental sides, the center value of Bd —> K + K~ (B s —> is a bit larger (smaller) than the data, 

which means there is a little room for us to search for possible effect of NP. Unfortunately, the CP asymmetries of 
these two decays have not been measured in the current experiments. 

3 The Contribution of The Z' Boson 

Now we shall study the possible contributions of the extra gauge boson Z' in these two decay modes. Ignoring the 
interference between Z and Z' bosons, we write the Lagrangian with Z' on the gauge interaction basis as 

C z ' = [MhPl + MhPr] ^ , (23) 

where the field ^ stands for the ith family fermion, g 2 for the coupling constant, e^ L (e^ B ) for the left-handed 
(right-handed) chiral coupling, and Pl,r = (1 T 7s)/2. When rotating to the physical basis, the mass eigenstates 
will be obtained by i/>l,r = V^l.r^l r j anc l the usual CKM matrix is given by Vckm = V UL Vj . Similarly, we can 
get the coupling matrices in the physical basis of up (down)-type quarks, 

= V ux e ux V,;l x , B* = V dx e dx Vl (X = L, R). (24) 

It is apparent that if £ u (d) L ^ R) are not proportional to the identity matrix, the nonzero off-diagonal elements in the 
B u ’ d will appear, which induce the FCNC interactions at the tree level. In the current work,we assume that the 
up-type coupling matrix tu L(R) are proportional to the unit matrix, and the right-handed couplings of are flavor- 
diagonal for simplicity. Thereby, the effective Hamiltonian mediated by the Z' , for example b —> sqq (q = u, d ) 
transition, is given by 

'B-eS = B bs(sb)v-A^2( B g q (qq)v- A + B f q (qq) v+A ) + h.c. , (25) 

where gi = e/(sin8w cos 8\y) and mz' is the mass of Z’ boson. The diagonal elements of the effective coupling 
matrices B BR are required to be real because of the hermiticity of the effective Hamiltonian. However, for the 
off-diagonal one of B z s , it might be a complex number and a new weak phase (f>b s is introduced, which might play 
important roles in explaining the large CP asymmetries in B —> Kn [19]. Since the above operators of the forms 
(sb)v-A(qq)v^A already exist in SM, we can represent the Z' effect by modifying the wilson coefficients of the 
corresponding operators. As a result, we reorganize the ea. (E5ll as 

n z s = -%VtbV: s E ( AC !+ a Cloi q) + A C^ + A CiOi q) ) + h.c.. (26) 
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Correspondingly, the contributions of the extra Z' boson to the SM wilson coefficients at the mw scale is given 

2 


AC ^ 3V tb V t * Vffi m z > 


ACo f7) =- 

3VtbV ts \gim Z ' 


f fJ‘ 2 m z \' 
\gim z >) 
f g2m z \' 
\gim z >) 


B, 


bs 


2 B L d { d R) ) , 


B, 


bs 


(B u { u R) 

(b^ R) - 5,^ 




(27) 

(28) 


One can see that the Z' contributes to the electro-weak penguins ACg ( 7 ) as well as the QCD penguins AC 3 ( 5 ). 
In order to show that the new physics is primarily manifest in the electro-weak penguins, we simply assume 


B 


L(R) 


~ —2 B dd R \ and this relation has been used widely uni munis]. Therefore, the Z' contributions to the 


wilson coefficients are 


where 


ac 3 s (5) = 0 , 


AC7| (7) = / ^tbV t * s \ ^ L { R ) e i4> bs j 


VtbV ts 


t LX = 


g2m z 

9\m Z ' 


jdL tdX 
^bs^dd 


VtbV t * s 


C X = L,R ), 


= Arg[S ( 


6s J 


(29) 


(30) 


Note that the other SM wilson coefficients at scale lower than m b will also receive contributions from the Z' boson 
through renormalization group (RG) evolution. Since in this model there is no new particle below mw, the RG 
evolution of the modified Wilson coefficients is exactly same as the one in SM [26] . 

Similarly, we also obtain the hamiltonian b —> dqq , and the corresponding wilson coefficients and inner functions 
are given as 


AG 3 d ( 5 ) = 0 , 


d _ /1 /-LL(R) ^it/> bd 


ACg (7) = 4 


(d X = 


( 92m z \ ‘ 

\gim Z ') 


tjL dX 

BbdB S s 


v tbV t * d 


td I M 

VtbV* d ^ 

C X = L,R ), 


<j) bd = Arg [Bfr d ] . 


(31) 

(32) 


LL LR 

Now, we are in a position to discuss the possible parameter spaces of C, s d ' and 4> bS:b d■ In particular, we 
assume 32 /ffi ~ 1 because we expect that both the hypercharge U{1)y gauge group and the extra 17(1)' have the 
same origin from some grand unified models. Furthermore, we also hope m z /m Z i is at the order of (^(lO” 1 ), so 
that the neutral Z’ boson could be detected at LHC experiment directly. Note that the mass of a leptophobic 
m Z ’ boson has not been constrained till now, as aforementioned in Sec[T] In addition, we need to determine 


the other parameters \B^ S \ 


\B^ d \, \B^ q \ and new weak phases 


1 bs with the accurate data from B factories and 


-rrO 


LHCb experiment. For example, B bd and (f> b s, bd could be extracted from B q -B q (q = d , s) mixing. In order to 


explain the mass differences between B q and with new Z’ boson, \B bs ^\ ~ \VtbV t *, d j\ is required. Then, with 
experimental data of B d , s nonleptonic charmless decays, B qq R ~ 1 could be extracted. For the new introduced 
phases 4> bs and <p bd , they have not been constrained totally, although many efforts have been done m, we therefore 
set them as free parameters. How to constraint of these parameters globally is beyond the scope of current work 
and can be found in many references muni ns]. So as to probe the new physics effect for maximum range, we 
assume £ ~ ( d R ~ Q d R £ [0.001, 0.02], ie, the range of m Z ' is about [636, 2800] GeV, and < j> bd ,bs £ [—180°, 180°]. 

In Fig(2[ we explore the possible effects of Z' boson on the decay mode B s —> ir + n~. In the left panel, we 
present the variation of the CP-averaged branching fraction as a function of the new weak phase <j> ba with different 
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Figure 2: The contribution of Z' to the decay mode B s —> 7r + 7r~. The left panel represents the branching 
fraction as functions of <j>b s , the dotdashed (green), dotted (red), and dashed (blue) lines represent results from the 
£ = 0.001, 0.01, 0.02, respectively. The region edged by dotdashed lines (black) is the experimental data, while 
edged by the solid lines (red) is prediction of SM. The right panel stands for the relation between the branching 
fraction and the direct CP asymmetry,the region edged by dotdashed lines (black) is the experimental data, while 
edged by the solid lines (red) is prediction of SM. 

£ = O.OOl(dotdashed), O.Ol(dotted), 0.02(dashed). The experimental region (filled by horizontal lines) and the SM 
predictions (filled by vertical lines) are also shown for the comparison. From this figure, one can see that the 
SM is consistent with the data within ler. Including the Z' contribution, it is apparent that the parameter space 
\(f)bs\ < 80° will be excluded. For |0& s | > 80°, if £ < 0.01, the contribution of Z' boson will be buried by the 
uncertainties of SM. One also sees that when £ = 0.02 the branching ratio will be enhanced to 7.6 x 10 -7 , which 
is larger than the SM prediction. Note that the averaged experimental have large errors, and the small band will 
help us to determine the magnitude of £. It is emphasized that the LHCb had obtained a bit larger result, which 
indicates the existence of a light Z' . In the right panel, we plot the relation between the branching ratio and the 
direct CP asymmetry A^jp with an extra Z' boson. The region edged by blue curve is the possible region with 
parameter £ < 0.02 and </>&„ € [—180°, 180°]. With the experimental data, the lower half of the region can be 
excluded. In the permitted region, the range of direct CP asymmetry is [—3.2%, 0.1%], which is much larger than 
the estimation of SM (the grey region). The future measurement of these values in LHCb (LHC-II) experiment 
and the high energy e + e _ collider will help us to probe the effects of Z'. Note that when discussing the effects of 
Z' boson, we will not include the uncertainties induced by wave functions and scale t. because the major objective 
of this work is searching for the possibility of new physics signal, rather than producing acute numerical results. 

Similarly, the effects of the extra Z' boson in Bd —> K + K~ are also presented in Fig. [3] In the left panel, it is 
clear that the position of SM prediction is on the top of the experimental data, although some parts of them overlap 
with each other. Furthermore, the heavy Z' contributions (£ < 0.1) are not apparent due to the uncertainties of 
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Figure 3: The contribution of Z' to the decay mode Bd —> K + K . The legends are the same as in Fig. [2] 

SM. Moreover, for the fod, the ranges [0°, 180°] and [—180°, —110°] will be excluded. We also plot the region (edged 
by curve) related to the direct CP asymmetry and the branching fraction, as shown in the right panel. Note that 
the SM prediction is 30% ~ 42%, but with a light Z' the estimated range is to be 5% ~ 60% after considering the 
constraint from the experimental data. It is concluded that for Bd —> K + K~ the future measurement of direct 
CP asymmetry will help us to search for the possible effect of a light Z', though its contribution to the branching 
fraction is polluted by the SM uncertainties. 

At last, we shall discuss the Z' effect on the CP symmetry parameters Sf and Tlf. For B s —> , as the 

weak phase of VtbV* a is very small, both Sf and "H/ are not sensitive to the NP. On the contrary, for Bd —> K + K ~, 
Sf and Hf are sensitive to the extra leptophobic Z’ boson. In Fig 01 we plot the relations of Sf (right panel) and 
Hf (left panel) with varying (f>bd from —180° to 180°, when ( = 0.01 and £ = 0.02. The estimations of SM (edged 
by the lines) are also presented. From the figures, one can see that in the permitted range of (f>bd, with a light Z' 
boson (£ = 0.02), Sf could reach —0.55, which is larger than the prediction of SM. For "H/, its values could reach 
to —0.75 when <pbd = —50°. The future measurement of them will help us to further constrain the parameters, 
which might helpful for direct searching for a light Z 1 boson. 


4 Summary 

In this work, we have studied the pure annihilation decays Bd —> K + K~ and B s —> 7r + 7r~ in the SM and the family 
non-universal leptophobic Z ' model. Although the SM predictions in the PQCD approach are in agreement with 
experimental data, the branching fraction of B s — > (Bd — > K + K~) is a little bigger (smaller) than the LHCb 

result, which may indicate the survival space for a light Z' boson. Inspired by this fact, we have constrained the 
U( 1) phase as 4>bd £ [—110°, 0°], and fos is |^b s | > 80°. Within the allowed space range, the direct CP asymmetry 
for B s —> 7r + 7r - is predicted as [—3.2%, 0.1%], while it is 5% ~ 60% for Bd —> K + K~. Furthermore, we have 
also calculated the mixing-induced CP asymmetries and found that the parameters Sf and %f of Bd — > K + K~ 
are very sensitive to the effect of Z’ . With a light Z\ the maximum (minimal) value of Sf(T~Lf) can reach —0.55 
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Figure 4: The CP symmetry parameters Sf (left panel) and Hf (right panel) as a function of the weak phase 0 m, 
the dotted (red) and dashed (blue) lines represent results from the £ = 0.01, 0.02, and the regions edged by solid 
line (green) are the predictions of SM. 

(—0.75). The future measurements of these observables may provide us some hints for direct searching for a light 
leptophobic Z' boson. 
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